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Based on results of a search for the lepton-family-number-violating decay K~* 



7r fi e 



with 



data collected by experiment E865 at the Alternating Gradient Synchrotron of Brookhaven National 
Laboratory (BNL), we place an upper limit on the branching ratio at 2.1 x lO^ 11 (90% C.L.). 
Combining the results with earlier E865 data and those of a previous experiment, E777, an upper 
limit on the branching ratio of 1.3 x 10 -11 (90% C.L.) is obtained. 



PACS numbers: 13.20.- 
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I. INTRODUCTION 



We report on a search for the decay K + — > 
7r + /i + e~ (if^jue). This is a lepton-family-number- 
violating (LFNV) decay and, thus, is strictly forbidden 
in the Standard Model (SM) with massless neutrinos. 
Incorporating massive neutrinos into the SM, which is 
required by the growing evidence for the lepton fam- 
ily transformations in the neutrino sector [1], results in 
a prediction of LFNV kaon decays at an unobservably 
low level [2]. Thus an observation of a LFNV process 
like the decay K niie would serve as a clear indication of 
physics beyond the SM. Moreover, extensions of the SM 
such as Extended Gauge theories [3, 4], Technicolor [5] and 
Supersymmetry[6] do allow LFNV processes. 

A first search for iVr^e, performed during the 1970s, 
resulted in an upper limit on the decay's branching ratio 
(B) of 4.8 x 10~ 9 (90% C.L.)[7]. A search for this decay 





FIG. 1: Left: Feynman diagram for Extended Technicolor 
mechanism for K + — » n + fi + e~ decay; right: Feynman dia- 
gram for K — > 7T° n + v decay. (0 C is the Cabibbo angle, gn 
and g are the ETC and weak coupling constants) 

in BNL was initiated by Experiment 777, in which the 
upper limit B < 2.1 x 10~ 10 was established [8]. The data 
collected in 1995[9] and 1996[10, 11] by BNL Experiment 



865, a successor to E777, allows us to lower the combined 
upper limit to B(K+ -> Tr+^+e") < 2.8 x 10~ n (90% 
C.L.) [11]. Nevertheless a null search result is usefull since 
it can either put constraints on the parameters of the 
existing extension models that allow LFNV or rule them 
out. 

For example, in the Extended Technicolor model 
(ETC) [3] the transition between the leptons of different 
generations can be mediated by a horizontal ETC boson 
(H). The corresponding diagram, shown in Fig. 1 is sim- 
ilar to the diagram of the familiar decay K + — > 7r°/i + iA 
Using the assumption that H and W boson couplings are 
approximately equal to each other, one can relate the H 
boson mass, Mu, to the K Viie branching ratio[3]: 



M H « 85 TeV 



10" 



B(K+ -> ?r+/z+e) 



1/4 



(1) 
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The quoted limit B(K+ -> ir+fi+e) < 2.8 x 10 
corresponds to a 65 TeV lower limit on the Mjj mass 
scale. A competitive constraint on the new physics mass 
scale may be obtained from another LFNV kaon decay, 
Kl —* [i^e^ (K^e), for which an upper limit on its 
branching ratio is set at 4.7 x 10~ 12 [12]. However, the 
fact that the decay is sensitive to the axial-vector 
and pseudoscalar (sd) quark transitions, while K^^ e is 
sensitive to vector, scalar, and tensor transitions, makes 
these two processes complementary in the search for new 
phenomena. 

This illustrates the significance of the rare decay 
searches which can be performed at an easily attainable 
low energy but at the same time can effectively probe the 
high mass region. 

The analysis outlined here is based on data recorded 
at the Brookhaven Alternating Gradient Synchrotron 
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FIG. 2: Overall view of the detector with a simulated K^^e event (horizontal cross section at beam height). C1,C2: gas Cerenkov 
counters; P1,P2,P3,P4: proportional chambers; D5,D6: dipole magnets; A,B,C,D - scintillation counter trigger hodoscopes. 



(AGS) during a five month run in 1998, employing the 
E865 detector. 



II. EXPERIMENTAL SETUP 



A. Overview and design philosophy 



To perform a successful search for a rare decay, like 
the K^fie, two basic conditions must be satisfied: first, 
the parent particle (K + ) must be produced in copious 
amounts and second, the detector must be able both to 
detect signal and to suppress backgrounds to a level low 
enough for a signal to be seen. The first condition was 
met by employing the high intensity proton source (AGS) 
to create the K + beam. With the projected K^^e sin- 
gle event sensitivity of 1CP 12 , in order to meet the second 
condition we designed the E865 detector to have excellent 
event reconstruction and particle identification (PID) ca- 
pabilities. The particular requirements for the detec- 
tor's event reconstruction and PID performance were de- 
rived from the careful consideration of the possible back- 
grounds. The dominant way to get a 7r + , /i + and from 
a kaon decay is through the decay chain K + — > 7r + 7r + 7r~ 
(K T ), tt + — > 7r~ — > e~V, which has a combined 

branching ratio of 6.8 x 10~ 6 . The fact that the in-flight 
pion decays have to happen within the detector and the 
neutrinos carry off momentum, distorting an event's kine- 
matics, led to the design of a high-resolution momentum 
analyzing spectrometer which allowed charged particle 
reconstruction and momentum measurement. The study 
of the Monte Carlo simulated events and off-line event 
reconstruction (to be discussed later) showed that such 



backgrounds could be suppressed to the level of 10 -13 by 
applying simple kinematic cuts. However, the K r decay 
would also mimic the K^^ e if a 7r~ was misidentified as 
e~ and tt + was either misidentified as [i + or underwent 
an in-flight decay (tt + — > H +v )- Another background 
was caused by rather common kaon decays: K + 
(K^), with a 7r + misidentified as and K + - 



7T° [l + V 



Both decays would mimic K^ae if the ir° un- 



derwent a Dalitz decay (ir° 



"7) and the e + was 



misidentified as a n + . In the following we will refer to 
both decays (K„2 and followed by a Dalitz decay of 
7T°) as Koai- The suppression of the K T and Koai back- 
grounds required a PID system with abilities to veto n~ 
(while still being efficient in identifying e _ ), veto e + and 
to discriminate fi + from it + . 

A detailed description of the apparatus (see Fig. 2) 
and the beam arrangement is presented elsewhere [13]. 
Here we give only a brief overview of the kaon beam and 
detector elements. 



B. Kaon beam 

The AGS at Brookhaven National Laboratory served 
as a source of the primary 25.5 GeV/c proton beam, and 
delivered about 1 x 10 13 protons to the kaon production 
target in a spill of 2.8 s with a cycle time of 5.1 s. Before 
the extraction from the AGS, the proton beam was de- 
bunched, i.e. particles were distributed over the circum- 
ference of the ring to obtain a uniform beam intensity 
during the spill. The extracted beam of protons, hit- 
ting a 15 cm long water-cooled copper target with trans- 
verse dimensions 5x5 mm 2 , produced a large variety 



3 



of particles, but mainly pions, protons and kaons sur- 
vive. After the production target, the unseparated sec- 
ondary beam was momentum selected at 6 GeV/c (±2%) 
and was transported to the beginning of the E865 detec- 
tor through the specially designed 28m long A2 beam 
line[13]. A major concern during the design of the A2 
beam line was the muon halo, produced by kaon and 
pion decays. It itly suppressed by tight colli- 

mation and two foci followed by dipole magnets, which 
effectively swept momentum-degraded muons out of the 
acceptance [13]. During data taking, the beam flux in the 
A2 beam line was estimated to be 2 x 10 8 K + , 4 x 10 9 
7T+, and 2 x 10 9 protons per spill. 



C. Detector overview 

The E865 detector was located at the end of the A2 
beam line, and its first element was the 5-m long evac- 
uated decay region within which about 6% of the en- 
tering kaons decayed. Downstream of the decay region, 
a dipole magnet (D5) swept the charged decay products 
away from the beam, with negative particles going mostly 
to the left, and positive to the right side of the appara- 
tus. The momentum-analyzing spectrometer consisted of 
proportional wire chambers (PWC), P1-P4 and a second 
dipole magnet (D6), which steered the particles back into 
the acceptance region of the detector elements located 
further downstream. The PWCs, each consisting of four 
wire planes, were desensitized in the region where the 
beam passed. This arrangement yielded a momentum 
resolution of <r p « 0.003 P 2 GeV/c, where P, the mo- 
mentum of the decay products in GeV/c, was typically 
in the range 0.6 to 3.5 GeV/c. 

The first part of the PID system consisted of two large 
atmospheric-pressure Cerenkov counters (CI and C2), lo- 
cated upstream and downstream of the P3 wire chamber. 
Each counter was separated by a thin membrane into 
two (left and right) parts providing independent particle 
identification for positive and negative tracks and allow- 
ing use of different gases. For the purpose of reducing 
the possibility of misidentifying a tt~ as e~ the left sides 
(C1L, C2L) were filled with a high Cerenkov threshold 
gas, hydrogen (74 = 60), and had a light yield of about 
2.3 photoelectrons (p.e.) for e~. To effectively register 
and veto e + , the right sides (C1R, C2R) were filled with 
a low-threshold gas, methane (j t — 35), and had a light 
yield of about 5.8 p.e. for e + . In order to reduce beam 
Cerenkov radiation in C1R and C2R, closed tubes filled 
with hydrogen gas were placed in the beam region. 

The next PID detector element was an electromagnetic 
calorimeter of the Shashlyk design[14]. It incorporated 
582 modules, 11.4cm by 11.4cm by 15 radiation lengths 
each, assembled in a 30x20 array with 18 modules in the 
middle removed for beam passage. Th e approxi mate en- 
ergy resolution for electrons was 8% / y/E(GeV). Typical 
energy deposition of a minimum ionizing particle was 250 
MeV. 



The last PID detector element was a muon detection 
system, that was located downstream of the calorimeter 
and consisted of 24 planes of proportional tubes inter- 
spaced by iron plates. The plate thickness was 5 cm be- 
tween the first eight pairs of planes and 10 cm between 
the last four. Finally, four arrays of scintillator counter 
hodoscopes, A,B,C and D, were used for triggering pur- 
poses. 



D. Trigger requirements 

The hardware trigger, described in detail fn[13], was 
designed as a four-level structure with increasing sophis- 
tication and response time at each level. The lowest trig- 
ger level (TO) selected events with three charged-particle 
tracks, one on the left and two on the right side, by re- 
quiring at least three coincidences between an A-counter 
slat and the calorimeter module behind it (A-SH). By us- 
ing a programmable matrix lookup unit (MLU), for each 
combination of coincidences between individual counters 
on the right, only a limited, kincmatically acceptable, re- 
gion on the left was allowed. The MLU was programmed 
to maximize the acceptance of the three charged body 
kaon decays by using simulated K n/le events. In order to 
reduce contamination from decays occurring downstream 
of the decay region, at least one coincidence on both left 
and right sides between the D-counter and A-SH was re- 
quired as well. The TO signal was ready in about 175 
ns after the particle hit the A-counter and had a rate of 
about 2.5 kHz. 

The next trigger level (Tl) used information from the 
Cerenkov counters and the muon range telescope for pur- 
pose of PID. In the case of K^^ e the dedicated trig- 
ger (MUE) demanded the presence of an electron and 
a muon. Consequently, the MUE trigger required sig- 
nals corresponding to at least 0.25 p.e. in both Cerenkov 
counters on the left side (C1L,C2L) to select e~ and a 
hit in the muon range telescope (B-hodoscope) to be spa- 
tially consistent with the A-SH hit in order to select /i + . 
The Tl signal was available about 130 ns after the TO, 
and had a rate of about 1.3 kHz. 

The next trigger level (T2) was designed to discrim- 
inate between events with high and low e~e + invari- 
ant mass. Information from coincident A-SH hits, and 
Cerenkov counters was used to select events consistent 
with low vertical separation of e + e~ tracks (low invari- 
ant mass). That trigger was not utilized in the K vlie 
search. 

The final trigger level (T3) compared the number of 
hits in each of the PWC plane with a predetermined like- 
lihood distribution obtained offline. This trigger rejected 
a small fraction of events that had either too few or too 
many PWC hits. During data taking, the average rate of 
T3 was about 1800 per 2.8 s spill or about 0.6 KHz. 

In addition to the MUE trigger, several prescaled mon- 
itor triggers were recorded, e.g. a minimum bias trigger 
(TAU), which required only TO signal and was dominated 
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by the K T events, and two triggers sensitive to the K^ai 
events. The first one (EEPS) demanded the presence of 
an e~ and e + by requiring signals (at least 0.25 p.e.) 
in both Cerenkov counters on the left and right sides 
(C1L,C2L,C1R,C2R) and was used for PID studies and 
background estimates. The second Koai sensitive trigger 
(CERENK) was designed to obtain an unbiased response 
of the Cerenkov counters and required only three out of 
four sides of the Cerenkov counters (C1L,C2L,C1R,C2R) 
to have signals above 0.25 p.e. 

During a five month data taking period in 1998 we 
recorded a total of 1.3xl0 9 triggers on tape, of which 0.3 
xlO 9 were MUE triggers, successfully reconstructed in 
the offline analysis. 



III. K nfle EVENT SELECTION AND OFFLINE 
ANALYSIS 

A. Event reconstruction 

The process of kinematic reconstruction of an event 
followed the TO hardware trigger requirements (A-SH) 
and started with the "clump" finding, i.e. determining 
positions of tracks in the calorimeter based on their en- 
ergy depositions and correlated A-counter hits. For each 
"clump" found in the calorimeter, a window in the PWC 
P4 was defined and space points, requiring at least three 
of the four wire planes in the chamber, were searched for. 
Likewise, the space points found in P4 determined win- 
dows in P3 and so forth. A track was formed if at least 
three PWCs contributed with a space point. Next, using 
the measured magnetic field map inside the dipole mag- 
nets, the momentum of each track was fit. For events 
containing at least three reconstructed tracks, a fitting 
algorithm was used to determine the decay vertex po- 
sition by minimizing the quantity S, the r.m.s. of the 
closest approach of the three charged tracks to the com- 
mon vertex point. For events containing more than three 
tracks the combination with the smallest S was tagged 
as the most probable set of track candidates from a kaon 
decay. Studies of fully reconstructed K T events showed 
that the mean value of S increased for vertices found 
further upstream (lower z) consistent with Monte Carlo 
simulations. The dependence could be fitted with a sec- 
ond order polynomial (denoted S(zj). The normalized 
quantity S norm — S/S(z), which is independent of z, 
was then used to judge the vertex quality in the K^^ e 
analysis. As the last step, the information from the PID 
detectors was assigned to each found track. 

B. Event Selection 

Basic requirements for K„^ e included the presence of 
a three-charged-track vertex within the decay region of 
an acceptable quality, S norrn . In addition, in order to 
ensure that reconstructed tracks came from a real kaon 
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FIG. 3: Target likelihood, Lrarget, distributions for the re- 
constructed K + — > 7r + 7r + 7r~ (upper) and K + — > 7r°7r + with 
a it" — > e + e~7 (lower) events. Poor target likelihood for 
the K + — ► tv°-k + events is caused by the photon, which is 
not included into the kinematic reconstruction of the kaon 
momentum vector at the vertex causing it to fall out of the 
acceptable phase space of the kaon beam. For signal events 
Lrarget > —20 was required, corresponding to an efficiency 
of 92%. 



decay with all the particles being detected, the recon- 
structed total momentum vector was tracked back to the 
production target using the first-order beam-line transfer 
matrix, and its phase space (the kaon momentum Pk+, 
the position x, y, and the direction, X , 6 y ) was com- 
pared to that of the kaon beam. The five-dimensional 
phase space of the beam was broken up into three two 
dimensional distributions: x versus 6 X , y versus 9 y and 
Pk+ vs x ■ The properties of the kaon beam were deter- 
mined from fully reconstructed K T events and were used 
as input for a corresponding logarithmic likelihood func- 
tion Ltarget (see Figures 3 and 9, and also Sec. HID). 
The invariant mass M wfle , calculated using the daughter 
particles masses {m n , m M , m e ) and their measured mo- 
menta, was required to be consistent with the kaon mass, 
493.67 MeV, within the resolution, a M ^ e ~ 4 MeV for 
simulated K^^ e events. 

Finally we required an unambiguous identification of 
e~, 7T+ and n + . Electron PID required a signal (at least 
0.3 p.e.) in the left sides of both Cerenkov counters 
(C1L,C2L) with corrected timing within ±5 ns. In addi- 
tion, the energy loss in the calorimeter (E) was required 
to be consistent with the momentum of the track (P), i.e. 
the ratio E/P to be between 0.65 and 1.5. Pions were 
identified requiring the absence of a signal above 1.2 p.e. 
on the right sides of the Cerenkov counters (C2R,C2R) 
and corresponding corrected timing of more than ±5 ns. 
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TABLE I: Identification efficiencies and probability of 
misidentification for the primary selection PID cuts. The 
symbol "— ►" stands for identified as. 



TABLE II: Main K + -> 7r + ^ + e _ backgrounds from other K + 
decays. 









— > e 


n+ 


0.804 ± 0.008 


0.066 ±0.001 










(8.7 ±2.6) x 10~ 6 






0.79 ±0.01 




e+ 


(1.1 ±0.1) x 10~ 5 


< (1.1) x 10~ 5 




e~ 






0.767 ± 0.003 



In addition, energy loss in the calorimeter was required 
to be consistent with that of a minimum ionizing particle 
or a hadron shower (E/P < 0.9). Muons were identified 
by requiring the absence of a signal above 1.2 p.e. in 
both Cerenkov counters (C1,C2), energy deposited in the 
calorimeter consistent with minimum ionization (E <450 
MeV) and penetration depth in the muon range consis- 
tent with its momentum. Table I summarizes the PID 
efficiencies and probabilities of misidentification achieved 
by the above mentioned selection cuts [15]. 

Although the selection criteria described above pro- 
vided a high sensitivity for the signal (respective effi- 
ciency of each cut was approximately 95%) they let a 
considerable number of background events pass. 



C. Overview of backgrounds 



Possible K nlie backgrounds could be classified into two 
different types. The first type consisted of rather com- 
mon kaon decays and was dominated by K T and Koai- 
These kaon decays could be tagged as K vlxe , if kine- 
matic requirements were satisfied and certain daughter 
particles were either misidentificd or underwent an in- 
flight decay as given in Table II. Due to the excellent 
particle identification capability of the E865 apparatus, 
this type of background could be efficiently suppressed 
by optimizing the PID selection in order to reduce re- 
spective misidentification probabilities, P(ir~ as e~) and 
P(e + as 7T+). The second type of background consisted 
of accidental combinations of the ir + , n + and e~ tracks, 
originating from separate kaon decays, that would sat- 
isfy kinematic requirements for K^^e- Such accidental 
events were usually characterized by uncorrelated track 
timing and poor kinematic quality (vertex quality S norm 
and Lxarget)- Instead of tightening the selection cuts to 
suppress these backgrounds, we chose to employ the like- 
lihood method, a brief overview of which is given below. 



K + Decay Mode 
K + — » 7r + 7r + 7r _ 
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K+ 

K+ 
with n 
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TrV^, 7T° 



e e 7 
¥ e + e~7 



None 



■n n as p T 

7r + as /i + and 7r~ as e" 
e + as 7r + and 7r + as /i 4 
e + as 7r + 



D. Overview of the Likelihood method 

In an analysis using cuts, all events that pass the fi- 
nal cuts are considered equally probable to be a signal. 
The likelihood method, however, allows one to test each 
event against a particular hypothesis, e.g. the likeliness 
for a particle to be of a certain type or the decay to be 
of a particular mode. For example a pion, passing all 
the electron PID cuts, would be misidentified as an elec- 
tron. However, for pions, the values for each of the PID 
variable will most likely fall at the edge of the various 
PID distributions while for a real electron they will be in 
the bulk. Comparing the probabilities of an electron and 
a pion having a particular PID response allows one to 
distinguish further between a real and misidentificd par- 
ticle. The same method can be applied to kinematic vari- 
ables. Thus, the use of multi-variable likelihood function 
allows one to differentiate quantitatively between signal 
and background on an event-by-event basis. 

The conventional definition of the log likelihood func- 
tion is 



£{x) = \ogP(x 



(2) 



where x — (xi,X2---x n ) is a vector of the measured kine- 
matic and PID response and P(x) is the n-dimensional 
probability density function (PDF). In the case where 
x\, X2---X n are independent, the probability can be sim- 
plified into: 

P(X!,X2, ...,x n ) = Pl( Xl ) ■ P 2 {x 2 ) ■ ... ■ P n {x n ) (3) 

where Pi(xi) is the ith one dimensional PDF, constructed 
from the corresponding distribution and x% is the respec- 
tive measured response. Using such an approximation, 
one can rewrite the log likelihood in the following way: 



£(f)=J>gPi(z0 



(4) 



i=l 
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FIG. 4: Left: M nfJ , e mass distribution for e + fj, + e~ events measured with the EEPS monitor trigger (Kn a l background). Vertical 
lines mark the M„- Me mass signal region; center: M ee invariant mass spectrum for measured Kr> a i events. The cut on M ee 
> 0.055 GeV removes 94% of the Koai events; right: M, F mass spectrum for the 7r + /i + 7r~ events measured with the TAU 
monitor trigger (K T background). Vertical lines mark the M^^ e mass signal region. 



Generally, an individual PDF can be represented as 
a continuous analytical function or a binned histogram. 
While it would be ideal to use the full range of a variable 
for the likelihood function to maximize the acceptance, 
it was more practical to use a finite (that would include 
around 95% of signal events) range, choosing particu- 
larly the most sensitive part of the distribution. In most 
of the cases the range of the likelihood function was de- 
fined by the values of the final cuts. The number of bins 
was selected to be large enough to reflect the shape of a 
particular distribution. Finally, in constructing a likeli- 
hood function, we required each PDF to be normalized 
to unity. 



IV. BACKGROUND STUDIES 

A. Kriai background 

A large part of the Kjj a i background was rejected by 
applying the invariant-mass cut (M ee > 55 MeV), since 
the invariant mass of the e + e~ pair from the K^ai decays 
is peaked at low values (see Fig. 4). Events containing 
photons registered in the calorimeter with a corrected 
timing within ± 2 ns, were rejected. Additional rejection 
was achieved by the Target Likelihood cut since, due to 
the presence of the photon that was not included in the 
kinematic reconstruction, the total kaon momentum vec- 
tor was usually inconsistent with that of the kaon beam. 
To estimate the level of the Koai background we used 
the EEPS triggered data. We selected e + fi + e~ events 
that could mimic K T ^ e (if an e + was misidentified as a 
7r + ) by applying the K vlie kinematic selection cuts and 
requiring identification of the /i + and e~ . Events, that 
were selected that way, were scaled by the EEPS trig- 
ger hardware prescale factor (20) and a corresponding 
misidentification probability - P(e + as 7r + ). To reduce 



the misidentification probability we constructed the log- 
likelihood function for the purpose of ir + identification. 
The likelihood included the ratio E/P and the quantity 
Emtio describing the spatial spread of the electromag- 
netic shower in the calorimeter. The latter variable is 
defined as the ratio of the energy deposited in the cen- 
tral module, i.e. the module to which a track points, 
to the total energy, which also includes the three most 
adjacent modules[15]. 

The selected cut value for the ir + PID likelihood, which 
resulted in the estimated misidentification probability 
P(e + as 7r + )=(2.5 ± 0.3) x 10 -6 , ensured a satisfactory 
suppression of the estimated K\]j a i background to the 
level of 0.09±0.01 events. 



B. K-r background 

Kinematically, the K T background was suppressed by 
applying the invariant M w ^ e mass cut, since assigning 
lower masses (rriejTO^) to the pions from the K T de- 
cay effectively removed mass from the system as seen 
in Fig. 4c. To estimate the level of the K T background 
we used the TAU triggered data. We selected n + [i + iT~ 
events that could mimic K^^e (if a ir~ was misidenti- 
fied as an e~) by applying the K^^e kinematic selec- 
tion cuts and requiring identification of the 7r + and n + . 
Events that were selected this way were scaled by the 
TAU trigger hardware prescale factor (10 4 ) and a corre- 
sponding misidentification probability - P(ir~ as e~). To 
reduce the misidentification probability we constructed 
the log-likelihood function for the purpose of e~ iden- 
tification. The likelihood included the number of p.e. 
and time registered by C1L and C2L, the ratio E/P in 
the calorimeter and again the variable variable E ra tio- 
The selected cut value for the PID likelihood, which 
resulted in the estimated misidentification probability 



7 




o n □ □ □□□□□□ □ □ □ n □ n r 

□□□□□□□□□□□□□□□□a 
□□□□□□□□□□□□□□□□a 




-10-8 



4 6 8 10 

K-K ("*) 



-v10 r 
1 8 
^6 

*** 4 ; 

2 

-2 
-4 

-6 : 

-8 
-10 



b) 



10-8 -6 -4 -2 



2 4 6 8 10 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

P*(GeV/c) 



*3 


0.16 


s= 






© 




0.12 




0.1 




0.08 




0.06 




0.04 




0.02 








d) 










n + /u,*e 












■' /•* ■ 











0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

P*(GeV/c) 



FIG. 5: Scatter plot of differences between track times: t w - 
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events [a)] and measured K T events[b)]; Scatter plot of ix + track angle at the vertex (ff) versus 7r + momentum (P n ) for the 
simulated K + — ► 7r () 7r + [c)] and K + — > 7r + /i + e~ [d)] events. 



P(7r as e )=(2.4 ± 1.4) x 10 6 , ensured a satisfactory 
suppression of the estimated K T background to the level 
of 0.1±0.1 events. 



C. Accidental background 

To reject accidental events effectively it was crucial to 
study and understand the precise accidental event mech- 
anisms. The simplest mechanism would be a combina- 
tion of two tracks from a real kaon decay and a single 
track from another decay, or beam background (2+1). 
Another possibility is that all three tracks would origi- 
nate from different decays or beam background (1+1+1). 
To determine what particular combinations of accidental 
tracks were dominant, the time structure of accidental 
events was examined. The time of each track (t e ,t w ,t^) 
was defined as the average of the times registered by the 
A counter and the calorimeter. The t e ~ t v and t e — 
track time differences were chosen to describe the event 
timing. A control sample of accidental three-track events 
was selected by requiring the reconstructed kaon momen- 
tum to be greater then that of the beam (P > 7 GeV/c) 
and removing the Lxarget cut. To increase the sample 
size, cuts on the vertex quality S norm , and the invariant 
mass M^e, were removed. Accidental events, selected 
in such a manner, showed a distinct cross in the scat- 
ter plot of the track time difference as seen in Fig. 5a, 
indicating that accidental events had two main compo- 
nents. The horizontal part of the cross identified the 7r + 
as the additional accidental track paired with the two 
tracks (e~/i + ) from another decay and the vertical part 
identified [i + as an accidental track. Accidental e~ tracks 
were not observed, since the scatter plot did not show a 
linear correlation between the variables. Finally, triple 
coincidences (1+1+1) appeared uniformly scattered and 
as we expected their contribution was negligible. For the 
signal events, the correlation between the track time dif- 
ference was examined using the measured response of the 
reconstructed K T events as can be seen in Fig. 5b. 

The difference in the time response between the signal 
and accidental events was used to create a time quality 



estimator variable, T max , which was defined as the max- 
imum absolute difference between the track times t n or 
tfi and the average times of the two remaining tracks t efl 
or i e7r respectively. The ratio between the two acciden- 
tal components was estimated to be approximately one 
to two (34% e^i+n and 66% e-K+^i). The most reason- 
able candidate for the two-track part (e7r) of the main 
accidental component (e7r + jj,) was the decay with 
a subsequent n° — > e + e~-f decay. In the two-body decay 
A%r2 the pion momentum and angle of the pion track at 
the vertex are correlated, as shown in Fig. 5c. On the 
other hand three-body decay signal events do not show 
such a correlation (Fig. 5d). The effect was best ob- 
served when the ir + angle-momentum dependence was 
projected on the axis orthogonal to the axis of almost 
linear correlation between the angle and momentum of a 
7T+ from the K^i decay as shown in Fig. 5c, equivalent to 
using the variable z = 6' w /0.14+P w /5 (P v in GeV/c and 
6 in radians) . By requiring the variable z to be less then 
one, 60% of accidental background events were rejected 
with 80% of them being from the en + n component. This 
caused only a 13.5% loss in signal sensitivity. 

The number of accidental background events was es- 
timated using the number of observed out-of-time K ntie 
candidates (3 ns < T max < 8 ns) which otherwise sat- 
isfied all K^fie cuts. To extrapolate into the acceptable 
in-time signal region, the timing distribution of the high 
momentum (P > 6.9 GeV/c) K^^e candidates was used. 
From that timing distribution we determined the pro- 
jection factor R (the ratio between the number of high 
momentum K v)le candidates with T max < 3 ns and the 
ones with the the T max in the range of 3 to 8 ns). The 
expected level of accidental background in the signal re- 
gion was thus determined to be 8.2 ± 1.9 events. Data 
study confirmed that the projection factor R, used for 
accidental background estimate, did not depend on the 
momentum cut-off value, which was varied in the range 
of 6.9 to 7.7 GeV/c. To test the validity of using high 
momentum accidentals to extract the time distribution 
of K vtie accidental background, we selected a sample of 
K-irne candidates with the M T(le mass cut reversed and 
no timing (T max ) cut applied. The timing response from 
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FIG. 6: The reconstructed M^^e mass distribution for the: a) estimated K T background; b) estimated Kr> a i background; c) 
estimated accidental background; d) Comparison of the M T(je mass distribution for the K nll& data (markers with error bars) 
and total estimated background (histogram). 



the latter sample was a superposition of the timing re- 
sponse from K^fie accidental background and the tim- 
ing response of the K T and JCdoZ backgrounds. The tail 
(T max > 3 ns) of the timing spectrum, which was domi- 
nated by accidental events in that sample, matched the 
timing distribution of the high momentum events. 



V. LIKELIHOOD ANALYSIS 



A. Check of the background estimates 

Following the principles of the so called blind analy- 
sis the Kir^e candidate events were examined, while the 
signal region was excluded by reversing the cut on the in- 
variant kaon mass {\M nfie - 0.4937| > 0.012 GeV). The 
invariant M I)ie mass distribution for the 216 observed 
K-Kfie candidates as illustrated in Fig. 6d exhibited a peak 
at lower mass values, which was caused by the K T back- 
ground. In order to check the validity of the background 
estimates, the M^^e mass distribution for the K v ^ e can- 
didates was compared with that of the predicted back- 
grounds. Following the same algorithm as that of the 
background estimates the M w/Je distributions were gen- 
erated for the main background processes using measured 
data: for the K T background, observed 7r~7r + /i + events 
were scaled by the misidentiheation probability P(jr~ 
as e~). For the Koai background, observed e~e + fi + 
events were scaled by the misidentiheation probability 
P(e + as 7T + ). Finally, for the accidental background, 
the out-of-time K ntie candidate events were scaled by 
the projection factor R (ratio of in-time to out-of-time 
accidental events). The satisfactory agreement between 
the predicted and observed events as displayed in Fig. 
6d served as a consistency check for our estimates of the 
background processes and misidentiheation probabilities 
(P(e + as 7T+) and P(n~ as e - )). 



B. General K vlle . likelihood 

The background study indicated that 98% of the back- 
ground to K vlxe was due to accidental events. As dis- 
cussed above, rather than imposing additional or tighter 
cuts, the likelihood function was chosen to suppress 
this background. The general K n/le likelihood was con- 
structed using the kinematic and timing variables which 
showed the most difference in their response for sig- 
nal and accidental events. The following six variables 
were selected: target likelihood, Lxarget, vertex quality, 
Snorm, reconstructed M I)ie mass, timing cut, T max , sum 
of the x 2 for the track reconstruction fits, and T extra , the 
minimum time difference between the time of the extra 
clump in the calorimeter and time of the track pair, which 
is closest in time. Data studies showed that within the 
statistical uncertainties the variables chosen were uncor- 
rected both for the signal and for the background. This 
allowed us to calculate the general likelihood for K^^ e 
and accidentals in the following manner: 

6 

Ci r p e {x 1 ...X 6 ) = ^log Pf^iXi) (5) 
i=l 

6 

C Acc ( Xl ...x 6 ) =^logPf cc ( Xi ), (6) 

i=l 

where X{ = {J-^Targett S norrn , Af^e, XjYacA:' ^maxi Textra) 

and P[^ e , pA cc are one dimensional PDFs for signal and 
background, respectively. For the signal, only the 
mass PDF was generated from the Monte Carlo simu- 
lated events with remaining PDFs generated from the 
measured data by using reconstructed K T events. For 
the accidental events, respective PDFs were generated 
from the data, using out-of-time K n/ie candidates, ex- 
cept for T max and T extra for which the PDFs were gen- 
erated using high momentum K vlxe candidates [15]. The 
total relative likelihood C for K n/J , e was dchned as the 
difference between C. Klie {xi...Xo) and Ca C c{x\---x§). 
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from the PDFs events for the K wfle (top) and the accidental 
background (bottom). The K^fie likelihood was calculated 
from the measured K T event sample (see text for details). 



The likelihood distribution for signal events was gen- 
erated using the measured data. With the exception of 
the invariant M niie mass, the response of the K T events 
for the likelihood variables was identical to that of the 
signal. The mass resolution of K T events had to 

be scaled by a factor of 1.725 in order to match the mass 
resolution of K nfle events. The latter was determined 
from the Monte-Carlo simulations[15]. Consequently, the 
signal likelihood template was generated using the good 
K T events with the scaled M n7T7T mass used as the M vlxe 
mass. 

The likelihood distribution for accidental events was 
obtained using a sample of out-of-time K viie candidate 
events. Since all six variables used in the likelihood had 
to be within the final cut values, the timing response 
was generated from a library of the measured timing re- 
sponse of the high momentum (P > 7 GeV/c) accidental 
events [15]. To test the assumption about the indepen- 
dence of the variables used in the general likelihood, their 
response was simulated directly from the corresponding 
PDFs. The fact that the simulated likelihood distribu- 



tions compared favorably with the ones derived from the 
data (see Fig. 7), confirmed our study of variables cor- 
relation and ensured the validity of the likelihood factor- 
ization of Eq. (3). 

C. Determining the number of K^^e events 
1. The Bayesian approach 

To determine the number of signal events, in the pres- 
ence of a significant number of accidental background 
events (8.2 ± 1.9), the Bayesian approach[16] was used. 
The likelihood function for Poisson-distributed data was 
defined as follows: 

L P (m;n) = \[ 1 (7) 

where 

k = arbitrarily chosen number of bins 
rii — the number of observed events in the ith bin 
n = (ni,n 2 , ...,nfc) 

fc 

N = total number of observed events (— ^ rii) 

i=l 

rrii = the number of expected events in the i t h bin 
to = (toi,to 2 , ...,TO fc ) 

The number of expected events in the ith bin of the 
likelihood distribution was parameterized in the following 
way: 

m,i = N sig P vfiei + N bg P acCi (8) 

where N sig was the number of K nlie events, N^g was the 
number of accidental events, and P^e; and P OCCi were 
the probabilities to be in the ith bin for the signal and 
the accidental background, respectively. The parameter 
Nb g was determined earlier from the study of accidental 
background. The probabilities (P-n^i and P acC i) were 
extracted using the likelihood templates, created for the 
K^ne and accidental hypothesis. The statistical uncer- 
tainty of the background estimate was accounted for by 
adding a Gaussian factor to Eq.(7): 

where 

= 8.2 and a bg = 1.9 

Under the parameterization of mj, given in Eq. (8), 
the Poisson-distributed data likelihood function, Eq. (9), 
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was transformed: Lp(m; ft) =>- L(N S i g , Nb g ). As the next 
step, the probability was normalized to unity: 



P(N sig ,N k 



L{N siq ,N bq ) 



J™dN bg f™L(N sig ,N bg )dN sig 



(10) 



To determine the upper limit on the number of signal 
events at the 90 % C.L. one integrates the probability 
P(N sig ,N bg ) from zero up to N 7 ™ x , until the 0.9 proba- 
bility level is reached: 



Ic r dN bgJo " 9 L(N slg ,N bg )dN s . 
Io° dN bg J °° L(N sig , N bg )dN. 



0.90 



(11) 



sig 



By solving Eq. (11) one can numerically determine the 
number of signal events at the 90% C.L. 



2. Signal region 

Eight events survived the final K vlxe selection, which 
was consistent with our background prediction of 8.2 ± 
1.9. Examination of the likelihood distribution of the 
surviving events demonstrates that they are clearly not 
consistent with the signal hypothesis (see Fig. 8). By 
solving Eq. (11) for the surviving events an upper limit 
on the number of K nlie events at a 90% C.L. was deter- 
mined to be N sig < 2.4. 

A simple cut on the likelihood would remove all sur- 
viving events with only a 3% loss in efficiency. 

The distributions of the six variables used in the gen- 
eral likelihood for the observed eight events along with 
the respective signal and background PDFs are presented 
in Fig. 9. 



VI. NORMALIZATION AND ACCEPTANCE 
CALCULATION 

To calculate the K^^e branching ratio, the total num- 
ber of K + decays had to be determined. The total 
number of K + decays can be determined indirectly by 
normalization to a decay mode that is well understood 
and observable in the detector, and has a relatively high 
branching ratio. We used the K T decay as the normaliza- 
tion mode. The upper limit on the branching ratio was 
determined according to the following formula: 



B(ir[ie) < B(Tnnr) 



N(ir[ie) Accep(iririr) 
N(7ttttt) Accep(TT/j,e) 



■C (12) 



where B denotes the branching ratio of the decay, 
N(n/j,e) is the 90% C.L. upper limit on the number of sig- 
nal events, N{-k-k-k) is the number of observed K T events, 
adjusted for prescale factors, Accep is acceptance of the 
detector system, and finally C is the correction factor ac- 
counting for efficiency differences between the selection 
of signal (K vlie ) and normalizer (K T ) decay modes (see 
Sec. VIC). 
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FIG. 8: Scatter plot of K n ^ e data (bold dots) and simulated 
signal(top) and background (bottom). The abscissa is the 
log-likelihood of the reconstructed events, the ordinate is the 
invariant mass (M n ^ e ) of the detected particles. The horizon- 
tal lines mark the 3<r mass region. 
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FIG. 9: Distributions of target likelihood (Lrarget), recon- 
structed kaon mass M^^e, vertex quality Snorm, minimum 
time difference between extra clumps in the calorimeter and 
closest in time track pair (T ex tra), sum of track's x 2 , arid T max 
for the Kn^e data. For the reference purpose corresponding 
PDFs are shown for signal (dotted line) and background (solid 
line) hypotheses. 
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TABLE III: Cuts used for determining the number of K T 
decays. 



Variable 



Description 



Efficiency 



Snorm 
3 

i=l 

LTarget 



iV 7 

S2trk 



vertex quality 

sum of track's fit \ 2 

target likelihood 

time quality indicator 

no extra clumps in the 
calorimeter within 1 ns 

photons veto 

two track vertex quality 

reconstructed kaon mass 



0.927 ±0.005 

0.967 ±0.005 

0.923 ± 0.005 
0.911 ±0.004 
0.956 ± 0.005 

0.965 ± 0.005 
0.959 ± 0.005 
0.976 ± 0.005 



A. Total number of K T events 

As mentioned in Sec. II D, we collected K T decays 
in a minimum bias TAU trigger concurrently with K vfle 
events. This allowed us to deduce the number of normal- 
ization K T events by analyzing the TAU triggered data. 
To reduce the systematic uncertainties in the acceptance 
ratio calculation, the normalization sample was kinemat- 
ically selected using cuts which were identical to the ones 
for K^^e- However, some of the cuts used for the K vfle 
selection, including the PID cuts, the invariant M ee mass 
cut, and a cut (z < 1) that rejected 7r+'s originating from 
Kn2 decays, were not applied. 

The application of the selection cuts (see Table III) 
yielded 77226 events. The total number of the K T events 
was determined after accounting for the hardware (10 4 ) 
and software (50) prescale factors: 

N(wwk) = 77226 x 50 x 10 4 = 3.86 x 10 10 events. (13) 



B. Monte Carlo simulations 

The simulation started with the kaon beam at the up- 
stream end of the decay region. The phase space of the 
beam was deduced from a large sample of K T events, 
for which the incident K + could be fully reconstructed. 
The K + was then allowed to decay in a preselected mode 
along its trajectory in the decay region. The interactions 
of the beam particles and the decay products with the 
detector were simulated using the GEANT3[17] software 
package from CERN. The simulated detector response 
was then analyzed using the same reconstruction and se- 
lection procedures as data events. 

In order to check the quality of the simulation, control 
distributions for various kinematic variables were com- 
pared to measured data for K T and K Dal decays. Due 
to the different phase space and different daughter parti- 
cles, those two decays have significantly different control 



distributions. The quality of the simulation is demon- 
strated in Fig. 10, which displays the vertex quality, 
Snorm, which is a crucial parameter for event reconstruc- 
tion; the position of the vertex along the direction of the 
beam, Z of the vertex, which depends on the detector ac- 
ceptance; and reconstructed kaon momentum, Pkaon, 
which is sensitive to the resolution. 

As an additional consistency check, we estimated the 
number of observed K T decays, as registered by the TAU 
trigger, from the number of observed low M ee mass K^ai 
events, as registered by the EEPS trigger, by apply- 
ing the acceptance ratio obtained from the Monte-Carlo 
simulations. For the calculation of the low M ee mass 
Koai events acceptance, we simulated K + — > 7r°7r + , 
K+ ttV + ^, K+ -> 7r°e+i/, K+ -> 7r+7r 7r decays 



(all with a subsequent decay of a 7r 



e + e 7) and 



weighted them by their measured branching ratios[16]. 
The predicted number of K T decays was consistent within 
the statistical uncertainty (1.5%) with the number of ob- 
served K T events [15]. This underlined our understanding 
of the geometrical acceptance and the efficiency of the 
various detector elements. 

For 2.5 x 10 5 simulated K T events 31600 were accepted, 
resulting in: 



Accep{inriT) = 0.1264 ± 0.0008 

10 5 sin 
cepted, resulting in: 



(14) 



For 2.5 x 10 5 simulated K^ ae events 17003 were ac- 



Accep{nne) = 0.0680 ± 0.0005 



and 



Accep(Tnnr) 
Accep(TT/j,e) 



1.86 ±0.02. 



(15) 



(16) 



To estimate the influence of alignment inaccuracies we 
varied the size of a dead (desensitized) region where the 
beam passed (see Sec. II C), which was the most sensitive 
parameter for both K nlie and K T decays acceptance. We 
considered three cases, for which we defined a smaller 
(Aa; = Ay =-0.5cm), nominal, and a larger dead region 
(Aa; = Ay =±0.5cm). While the changes in the accep- 
tance of the K T and K vtie decays were noticeable, the 
acceptance ratio showed the following variation: 

• 1.88 ± 0.02 - smaller dead region; 

• 1.86 ± 0.02 - nominal dead region; 

• 1.83 ± 0.02 - larger dead region. 

From this study we deduced a systematic uncertainty 
of roughly 1.5% on the acceptance ratio: 



Accep(irinr) 
Accep(ir[ie) 



= 1.86 ± 0.02(stat.) ± 0.03(syst.) (17) 
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FIG. 10: Comparison of the Monte Carlo simulations of K T decays (histogram) with data (markers with error bars). Left: 
distance of closest approach S n0 rm to the common vertex for the three charged tracks; center: distribution of decay vertices 
along the beam direction Z (Z = at the entrance of the first dipole magnet); right: Reconstructed K + momentum. 



TABLE IV: Efficiencies used to calculate correction factor. 

Effic. Description Value 

1~ 7T+ detection 0.781 ± 0.004 

e M n+ detection 0.79 ± 0.01 

e e e" detection 0.745 ± 0.003 

eM ee invariant mass cut M e - V + > 55 MeV 

and M e - M+ > 55 MeV 0.919 ± 0.007 

€cl e~ required on the left side 

of the Cerenkov counters 0.926 ± 0.008 

ecR f + required on the right side 

of the Cerenkov counters 0.891 ± 0.008 

e z reject tt+ from K + -> 7r°7r + decay 0.865 ± 0.009 



C. Correction factors 

The correction factor C in Eq. (12) contains efficien- 
cies to account for the difference in reconstructing K T 
and -ftT^e events. As described in section VIA, most 
of the kinematic cuts for the K T sample selection were 
also applied to select events and the corresponding 

efficiencies were the same for both decay modes. The cor- 
rection factor was calculated, as a reciprocal of the the 
product of the selection cut efficiencies specific to the 
K niie (see Table IV), to be C = 3.32 ± 0.08. The bulk 
of C, 2.17, is the reciprocal of the product of the 7r+, 
/i + and e~ PID efficiencies, while the remaining factor of 
1.53 results from the specific K n/Jie kinematic cuts. 



VII. RESULTS AND CONCLUSIONS 

With the B(iririr) = 5.58 ±0.03% [16], and all the other 
factors in Eq. (12) determined, the upper limit on the 
K^ne branching ratio is set at: 

B(K+ -► TT+fi+e-) < 2.1 x 10- 11 (90% C.L.) (18) 

The details of the analysis of the final data set ob- 
tained by E865 are described in this paper. Here, we 
summarize the search for the decay K + — > 7r + /i + e~ at 
BNL, which included a series of four measurements in 
which the following 90% C.L. upper limits were obtained: 
2.1 x 10- 10 [8], 2.0 x 10- 10 [9], 3.9 x 10- n [10, 11], and 
2.1 x 10~ n (this paper). Since the last two measure- 
ments were not background free, we analyzed the com- 
bined likelihood function constructed as the product of 
the likelihood functions for the separate measurements 
to establish a combined upper limit. For the final data 
set, this function is described above in Eq. 9. The anal- 
ogous function, based on the distributions shown in Fig. 
4 of Ref. [11] was used for the 1996 data. The only com- 
mon parameter was a K + — > 7r + /i + e~ branching ratio B, 
introduced into the likelihood function via 

N ■ — N e ff B 

where N^* was the effective (corrected for K^^e accep- 
tance) number of kaon decays in a given experimental 
run. For the background and signal free measurements 
of Refs. [8] (E777) and [9], the likelihood function (Eq. 
9) is reduced to the simple form exp (—N^/^ B). 

Using the Bayesian approach in the analysis of the 
combined likelihood function, we have obtained 

B(K+ tt+m+O < 1.3 x 10~ n (90% C.L.) 

Although no evidence of new physics was found, the 
parameter space of the existing extension theories that 
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allow LFNV was reduced. Particularly, as discussed in 
Sec. I, the mass, Mjj, of the corresponding Extended 
Technicolor boson, with strength equal to that of the 
weak interaction, is bound by the limit obtained on the 
K^fie branching ratio to be greater then w 80 TeV. 
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